Zero-bias anomaly in disordered wires 
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We calculate the low-energy tunneling density of states u[e,T) of an A''-cliannel disordered wire, 
taking into account the electron-electron interaction non-perturbatively. The finite scattering rate 
l/r results in a crossover from the Luttinger liquid behavior at higher energies, v cc e'^ , to the 
exponential dependence i'{e,T = 0) oc exp(— e*/e) at low energies, where e* oc 1/(Nt). At finite 
temperature T, the tunneling density of states depends o n the energy through the dimensionless 
variable e/^/e*T. At the Fermi level v{e = 0,T) cx exp{-^e* /T). 
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The influence of electron-electron interactions on 
transport in disordered systems has been extensively 
investigated for the past two decades It is well 

known that the interaction has the strongest effect in 
low-dimensional systems. Electron tunneling into a one- 
dimesional conductor is suppressed by interactions even 
in the absense of disorder. This suppression, which 
yields vanishing tunneling density of states (TDOS) at 
the Fermi level, can be described in the framework of the 
Luttinger liquid theory. The recently discovered carbon 
nanotubes provide a unique opportunity for studying in- 
teraction effects in quantum wires Although the 
properties of single-wall nanotubes are well described by 
the Luttinger liquid theory the transport in multi- 
wall nanotubes (MWNT) is still not very well under- 
stood. The number of channels, disorder strength, and 
carrier concentrations in these systems can vary over a 
wide range and are difficult to control experimentally. 
Whereas some measurements indicate ballistic electron 
transport M , most experiments exhibit diffusive electron 
motion p|,p[-JlO[. Furthermore, experiments jll] demon- 
strate a strong suppression of TDOS v{e) near the Fermi 
level (e = 0). On the other hand, the existing microscopic 
theory treats the screened Coulomb interaction in the 
first order of perturbation theory; it provides the result 
for the correction to the density of states, 5v{e) oc —Xj^/t, 
which is valid as long as Sv[t) is small. Clearly, this result 
of the lowest-order perturbation theory is insufficient for 
finding the behavior of TDOS v(t) in the limit of e ^ 0. 

In this paper we present a theory of the zero-bias 
anomaly in the tunneling density of states in quantum 
wires. We treat the dynamically screened Coulomb in- 
teraction non-perturbatively and allow for an arbitrary 
value of er, where r is the elastic momentum relaxation 
time of electrons, and energy e is measured from the 
Fermi level. This enables us to describe the crossover 
from the known Luttinger liquid results ||l2| valid at 
higher energies, er 3> 1/\/N, to the new low-energy 
(er < I/VN) behavior of the TDOS: 



v{e, T) cx exp 
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The 



(hereinafter we use units with h — ks = 1) 
characteristic energy e* here depends on the interaction 
strength g, 



g Tre^ d 
, q — In — 



(2) 



on the number of channels N in the quantum wire, and 
on T. Here v is the Fermi velocity averaged over all chan- 
nels, and d 3> i? is the distance at which the electric field 
is shielded [p^ say, by conducting electrodes surrounding 
the wire. The scaling function F{x) and its asymptotics 
at X ^ 1 (regime considered first in Ref. 10) and x:$> 1, 
are presented in Fig. 1. The result (|l|) applies at suffi- 
ciently low temperatures and energies, when the value of 
the exponent in this equation is large. Note that accord- 
ing to Eq. (|l]), at finite T the characteristic scale for the 
energy-dependence of TDOS is given not by T, but by a 
much larger value %/ e*T. 




FIG. 1. The scaling function F{x) and its asymptotics: 
F{x) = 1.07 - x/2 for a; < 1 (dotted line), and F{x) ~ 1/a; 
for a; 2> 1 (dashed line). 
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In the absence of shielding Eq. is somewhat modi- 
fied, as the interaction parameter g becomes a weak func- 
tion of temperature and energy. We briefly discuss this 
case in the end of the paper, see Eqs. (p^)-(^l|). 

Having in mind experimental applications, we focus 
our discussion on carbon nanotubes. A typical MWNT 
consists of several (up to ten) graphite monolayer sheets 
rolled concentrically into cylinders. At zero doping they 
can be either metallic or semiconducting, depending on 
the helical arrangement of the carbon hexagons. The 
electron band structure of a single carbon nanotube p^ ] 
has two points in the Brillouin zone with the Dirac- 
like spectrum efe — v{k'^ + fc^)^/^; the velocity here 
is w ~ 8 • 10^ cm/s, and the transverse momentum 
k± = n/R is quantized due to the periodic boundary 
conditions around the circumference, —kpR < n < kpR. 
The number of conducting sub-bands around each Dirac 
point N = 2k pR is determined by the radius of the nan- 
otube and by the doping level, /i — vkp (~ 0.5 eV). 
A typical radius of the outermost shell is of the order 
of 10 nm. Each sub-band has its own Fermi velocity 
Vn — vy^l — n'^ /kpR'^ and momentum kpVn/v along the 
cylinder axis. The electrons are scattered between differ- 
ent sub-bands within the same tube by impurities, lattice 
imperfections and by the incommensurate lattice poten- 
tial of the neighbouring tubes. We consider the exper- 
imentally relevant case of the ballistic electron motion 
around the circumference, I = vt > R, and concentrate 
on low-energy, e < v/R, limit. Electron tunneling at 
higher energies, e > vl/R"^, in the opposite case of a 
short mean free path, I < R, was recently discussed in 
Ref. @. 

In the measurements of the TDOS the tunneling cur- 
rent propagates through the outermost shell ^ while the 
inter-shell tunneling is largely suppressed. The role of the 
electrons in inner shells is then believed to be reduced 
merely to the dynamical screening of the Coulomb in- 
teraction between the electrons in the outer shell. Little 
is known about the contribution of inner-shell electrons 
to screening since the doping level in the inner shells is 
difficult to characterize experimentally. Two distinct sce- 
narios can be imagined: i) the dopants are outside the 
nanotube, and the doping electrons reside in the outer 
shell only, ii) the dopants are distributed uniformly in- 
side the MWNT, which leads to a uniform density of 
carriers across the shells. Below we concentrate on the 
first scenario, when the inner shells of the tube may be ig- 
nored. The second scenario will be considered elsewhere. 

We start with calculating the TDOS in the first or- 
der in the screened interaction potential. This calcu- 
lation follows the well-known route first developed for 
the case of a diffusive electron motion |17|, and extended 
later to the case of an arbitrary value of er. The 

zero-temperature result can be cast in the familiar [ pT|j20| 
form. 
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(3a) 
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Here F is the impurity-renormalized vertex. Its inverse 
is given by the usual impurity ladder, 



F-i(c.,q)^ 



uj — qv 



■ LU — qv + i/r 



}_Sr^ _i w(w + i/r) - ((ju„ + q_LW_L„)2 

7yn ^ ^ 



(4) 



The product of Green functions averaged over the Fermi 
surface [abbreviated as in Eq. (|3a|)] equals 



e'(^,q) 



1 



[ui — qv + i/tY ^ 



where qv = qVn + qA_vj_n and v±n = vn/{kpR) is the 
transverse velocity in the n-th band, vq — X]ri('''^")~^ 
is the total density of states in the outermost shell (the 
summation accounts also for both spin directions and the 
presence of two Dirac points in the Brillouin zone). 

The function U{uj,q) in Eq. (3b) represents the dy- 
namically screened Coulomb interaction of electrons and 
is given by 



U{co,q) = 



^(q) 



l-V^(q) n(cc>,q)- 



(6) 



Within the assumptions of our model, the polarization 
operator here. 



n{uj,q,qm) = iyoT{uj,q) 



qv 



'uj — qv + i/r 
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is provided by the outer-shell electrons. The bare 
Coulomb potential in Eq. m) is 



V{q, q„ 



2e2 



d(f)Kt:, ( 2gi?sin— I cosmc/), (8) 



where Kq{x) is the modified Bessel function, and we used 
the fact that the momenta along the circumference of the 
tube are quantized and given by qm = m/ R. 

To consider the low-energy behavior of Si'^e), we will 
need only the long-range limit, gi? ^ 1, of Eqs. (|1)-(||); 
we find for the bare interaction 



V{q, qn 



- ) l™l' 



e^ln [min(|7, ^)], 



TO 7^ 0, 
TO = 0. 



(9) 



The integrals in Eqs. ( |3a| ) and (^) are dominated by the 
region of high frequencies cu ^ qv, where the dynamically 
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screened Coulomb interaction (|^) has plasmon poles, and 
can be rewritten as 



dq {uj + i/T)V{q,q.m) 
-oo 27r2 u [uj{uj + i/r) - ^^(q)] ' 



(10) 



Here ijJm{q) denotes the frequency of the plasmon exci- 
tations which in the long wave length limit, qd <^ 1, is 
given by the following equation: 



e^i^o^^wfl ln[min(-^, ^^i^)], m = 0, 

TO 7^ 0. 



|m| \1 ^11 + _R2 



(11) 



In Eq. we introduced the average squares of the lon- 
gitudinal, V?,, and transverse, electron velocities 



J2n 



For e > vj_/R in Eq. (3a), we can approximate the 
sum over to in Eq. (M) by an integral recovering the 
ballistic counterpart |19(| of the two-dimensional diffu- 
sive correction discussed by Egger and Gogolin jl^ for 
short-range interaction. However, in contrast to their 
conclusions, for lower energies e < v±/R the contribu- 
tion of TO ^ terms becomes energy independent for 
both short-range and Coulomb interaction, in the latter 
case because of the gaps in plasmon spectra. The m = 
plasmon, on the other hand, is gapless. Its contribution 
to Eq. (|o|) depends on e, and has a singularity at e ^ 0. 
Therefore, to study the energy dependence of DOS at 
low energies we neglect the non-singular contribution of 
the rnj^ modes and retain only the m — term in 
Eq. (|lo|). The expression ( ITTI) for m — ceases to be 
correct at frequencies larger that \/Nv/R which corre- 
spond to plasmons with wavelength of order of the tube 
radius R, representing the obvious ultraviolet cut-off for 
one-dimensional effects. Performing the integral over the 
momenta in Eq. ([lo|), we obtain 



V{io) = - 
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(12) 



where v ~ "Ylm ^n/N is the average Fermi velocity. The 
two distinct regions of the frequency dependence here, 
uj ^ 1/r and w <C l/r, define two domains for the 
energy-dependent correction to the TDOS. Substituting 
V{uj) into Eq. (|3a|), and assuming that vt > d/^/N, we 
obtain 



2^0 




25 ^^+lnAr, 
N 1 ln(A/e), 



e < l/r, 
e > l/r. 



(13) 



where g is defined in Eq. (|), and A = v\fN / [R^d^/^ . 
The first term of the low-energy asymptotic here is famil- 
iar from the zero-bias anomaly theory of Altshuler and 



Aronov jl] , while the second term represents the omitted 
in ||l[ high-frequency {ui ^ I/t) contribution to the inte- 
gral Eq. (3a). The behavior of (5;^ at e ^ 1/t corresponds 
to the ballistic electron motion. 

The perturbative expressions diverge at the Fermi level 
(e — > 0). To calculate the tunneling DOS to all orders 
in the interaction constant at small e, we make use of 
the phase approximation for the fluctuatingpotential in- 
duced by the electron-electron interaction Q . The non- 
perturbative expression for the density of states can be 
cast in the form equivalent to the one derived in Ref. pi| , 



^ ' ' ^ T cosh — 
2T 



X exp 



dt 



cos et 



^ ^ cosh^ 
sinh 2y 



cosh TiTt 
- cos Lot 



(14) 



In the region of validity of Eq. (|l|), the time integral can 
be evaluated within the saddle-point approximation. As 
the saddle point lies on the imaginary axis in the interval 
< —it < 1/2T, the denominator in Eq. (14) is always 
a slowly varying function and does not contribute to the 
saddle point exponent. 

We consider first the ballistic regime, e > ^ t* Jt. It 
allows us to reproduce the conventional Luttinger liquid 
results, so for brevity we mention here only the limit 
T ^ 0. Evaluating the integrand of Eq. dlj), we find 



a ■ 



(15) 



In the diffusive regime, e < ye%^, the main contri- 
bution to the TDOS comes from the small frequencies, 
where one can reduce Eq. (|l|) to V(w) = 7^*7 

TTLO^ . Fur- 
thermore, if e ^ 1/Nt and T ^ g/irNr, the saddle 
point approximation is again applicable for the evalua- 
tion of the time integral in Eq. (|l^. One can easily 
check that at the saddle point the coshTrTt function in 
Eq. (|l^ ) can be replaced by 1, and therefore the density 
of states satisfies the scaling form, 



K6,T) 



exp 



(Ar)" 
Here the function 



coshf - cosh?/Z5(x) 
F{x) = / dy ^ 




y3/2 gjjlJl I 



(16) 



x[i-z.(x)] (17) 



is determined by the value of the integrand in Eq. (|l4| ) 
at the saddle point tg = iZs/T . The dependence of Zs on 
the ratio e/V €*T = x is given by equation: 



1 dy sinhzsj/ 



(18) 
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Numerical solution of this parameter-free equation and 
the subsequent evaluation of the integral in Eq. ( |l7| ) 
yields the graph of scaling function F{x) plotted in Fig. 1 
and its asymptotics given in the figure caption. Note that 
the pre-exponential factor in Eq. ( p^ , which we omitted 
in Eq. (|l|) for the sake of brevity, provides the proper 
matching of the resu lts (p^ ) and ( p6| ) obtained in the en- 
ergy domains e ^ a/e* /r and e ^ ■\/ e* /t, respectively. 

We derived our main results, Eqs. assum- 
ing that the Coulomb interaction is shielded at distances 
~ d. In the absence of shielding, the long-range nature 
of the interaction potential leads to a stronger |Q than 
predicted by Eq. ( |l5|) suppression of the TDOS in the 
ballistic regime, ln[i^(e)/t'o] oc — [ln(e/e*)]'^/^. At lower 
energies (corresponding to the diffusive regime), the ef- 
fect of the long-range potential can be accounted for by 
replacing the parameter e* of Eq. (||) with the logarithmic 
function of energy and temperature. 



e*(e,T) -eoln- 



v/R 



4:VNT 



(19) 



After the definition of e* is adjusted to reflect this re- 
placement, we can use Eqs. ([r^)-(|l^) again. In the limit 
of low temperatures, T <C e^/eo, we find: 



i^(e) (X exp In — 

^ e eR 



(20) 



The suppression of the TDOS near the Fermi surface 
(e ^ %/Teo) at finite temperatures is given by 



for the TDOS measurement should be attached to the 
MWNT within a distance shorter than the localization 
length from each other. 

We presented above the derivation of TDOS only for 
the case of a relatively long mean free path /, exceeding 
the radius i? of a MWNT, but our main result (|^) is valid 
for any relation between / and R. The only additional 
feature appearing in the case Z <C i?, is the studied in 
Ref. ||l^ intermediate range of energies, where the TDOS 
behaves as in a two-dimensional disordered conductor. 

In summary, we have obtained the tunneling density 
of states in a disordered quasi-one-dimensional conduc- 
tor. Our results are non-perturative in the electron- 
electron interaction, and cover both the diffusive and 
ballistic regimes of the electron motion. In contrast to 
the two-dimensional case [21 23, the non-perturbative 



results ( 14 ) and ( pO[ ) are not given by a simple exponen- 
tiation of the first-order interaction correction (^3|). 
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The origin of the strong sup pression of the TDOS, see 
Eqs. (|l|), (111), (lol), and (|l^, lies in the re-distribution 
of charge of the tunneling electron along the wire. This 
process is impeded by the finite propagation time of 
plasmons which enable the charge spreading. The re- 
quirement that the wavelength of the relevant plasmons 
v^/N /e be shorter that the length of the wire L imposes 
the lower energy limit for the applicability of the present 
theory. For smaller energies, e < v\/N/L, the tunnel- 
ing DOS of the wire depends on the impedance of the 
leads attached to the segment To observe a sizable 
suppression of TDOS, the nanotube must be sufRciently 
long, and therefore have high intrinsic resistance. We find 
that in order to reach the strong diffusive renormalization 
of TDOS anomaly (e, T <C e*) the total resistance of the 
segment should be made larger than (h j e^) \J~Njg . This 
condition on the overall segment length does not inval- 
idate the employed method which assumes the diffusive 
motion of electrons. Indeed, the characteristic frequen- 
cies of the plasmons involved are high enough for the 
weak localization's corrections to be ignored in the cal- 
culation of TDOS. In order to avoid localization effects in 
the DC transport measurement, the two junctions needed 
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